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Abstract Determining the coronal electron density by the inversion of white-light polar- 
ized brightness (pB) measurements by coronagraphs is a classic problem in solar physics. 
An inversion technique based on the spherically symmetric geometry (spherically symmet- 
ric inversion, SSI) was developed in the 1950s and has been widely applied to interpret 
various observations. However, to date there is no study of the uncertainty estimation of this 
method. We here present the detailed assessment of this method using a three-dimensional 
(3D) electron density in the corona from 1.5 to 4 R Q as a model, which is reconstructed 
by a tomography method from STEREO/COR 1 observations during the solar minimum in 
February 2008 (Carrington Rotation, CR 2066). We first show in theory and observation that 
the spherically symmetric polynomial approximation (SSPA) method and the Van de Hulst 
inversion technique are equivalent. Then we assess the SSPA method using synthesized pB 
images from the 3D density model, and find that the SSPA density values are close to the 
model inputs for the streamer core near the plane of the sky (POS) with differences generally 
smaller than about a factor of two; the former has the lower peak but extends more in both 
longitudinal and latitudinal directions than the latter. We estimate that the SSPA method may 
resolve the coronal density structure near the POS with angular resolution in longitude of 
about 50°. Our results confirm the suggestion that the SSI method is applicable to the solar 
minimum streamer (belt), as stated in some previous studies. In addition, we demonstrate 
that the SSPA method can be used to reconstruct the 3D coronal density, roughly in agree- 
ment with the reconstruction by tomography for a period of low solar activity (CR 2066). 
We suggest that the SSI method is complementary to the 3D tomographic technique in some 
cases, given that the development of the latter is still an ongoing research effort. 
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1. Introduction 

The electron density is a fundamental parameter in plasma physics. Knowledge of the 
three-dimensional (3D) electron density structure is very important for our understand- 
ing of physical processes in the solar corona, such as coronal heating and acceleration of 
the solar wind (e.g., Munro and Jackson, 1977; Cranmer et al., 1999). The density struc- 
ture of the corona strongly affects the propagation of CMEs (Odstrcil and Pizzo, 1999; 
Riley, Linker, and Mikic, 2001; Odstrcil et al., 2002; Manchester et al., 2004). The den- 
sity is also important for estimates of the Alfven Mach number and compression rate of 
CME-driven shocks (Reames, 1999; Sokolov et al, 2004; Manchester et al., 2005), and 
for the interpretation of solar radio emission such as type II and type IV radio bursts pro- 
duced by coronal eruptions (Caroubalos et al., 2004; Cho et al, 2007; Shen et al., 2013; 
Ramesh et al., 2013; Zucca et al., 2014). 

The K corona arises from Thomson scattering of photospheric white light from free elec- 
trons (e.g., Billings, 1966). Because the emission is optically thin, the measured signal is 
a contribution from electrons all along the line of sight (LOS). The derivation of the elec- 
tron density in the K corona from the total brightness ( B ) or polarized radiance ( pB ) is a 
classical problem of coronal physics, first addressed by Minnaert (1930) and Van de Hulst 
(1950). Because of difficulties in generally separating the K-coronal component from the F- 
coronal component arising from interplanetary dust scattering (e.g., Billings, 1966; Hayes, 
Vourlidas, and Howard, 2001), most of the inversion techniques are useful for pB measure- 
ments. The F-coronal polarization is not very well understood (see reviews by Koutchmy and 
Lamy, 1985; Kimura and Mann, 1998). It is generally accepted that the polarized contribu- 
tion of the F corona can be ignored within 5 R e (Mann, 1992; Koutchmy and Lamy, 1985; 
Hayes, Vourlidas, and Howard, 2001), but some observations show that the F corona is 
almost unpolarized even at elongations ranging from 10 to 16 Rq (Blackwell and Pet- 
ford, 1966a, 1966b; Blackwell, Dewhirst, and Ingham, 1967). In addition, the F corona 
dominates the total brightness of the corona beyond about 4 or 5 R Q and makes it dif- 
ficult to recover the much fainter K-corona emission (Saito, Poland, and Munro, 1977; 
Koutchmy and Lamy, 1985; Hayes, Vourlidas, and Howard, 2001). To retrieve the elec- 
tron density of the corona from a single 2D pB- image, one needs to assume some spe- 
cial geometries for the distribution of electrons along the LOS. The previous studies have 
modeled the electron density distribution in several ways, including the simple spher- 
ically symmetric model (Van de Hulst, 1950), the axisymmetric model (Saito, 1970; 
Munro and Jackson, 1977; Quemerais and Lamy, 2002), or the models that take into ac- 
count large-scale structures, such as polar plumes in the coronal holes (Barbey et al., 2008), 
or active streamers in the equatorial regions (Guhathakurta, Holzer, and MacQueen, 1996). 

Among the above inversion methods, the spherically symmetric inversion (SSI) devel- 
oped by Van de Hulst (1950) (called the Van de Hulst inversion thereafter), was the most rep- 
resentative and commonly used. He found that the density integral for pB signals becomes 
invertible if the latitudinal and azimuthal gradients in electron density are weaker than the 
radial gradient (i.e., a local spherical symmetry approximation). This classic inversion tech- 
nique has been applied to establish the standard density models of the coronal background 
at equator and pole in the solar minimum and maximum (see Allen, 2000) and the density 
models of near- symmetric coronal structures such as streamers and coronal holes (Saito, 
Poland, and Munro, 1977; Gibson and Bagenal, 1995; Guhathakurta and Fisher, 1995; 
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Guhathakurta, Holzer, and MacQueen, 1996; Gibson et al., 1999). The SSI method was 
also used to analyze detailed density distribution of fine coronal structures observed in 
eclipses ( e.g ., Koutchmy, 1994; November and Koutchmy, 1996), and to derive the 2D den- 
sity distribution of the entire corona (Hayes, Vourlidas, and Howard, 2001; Quemerais and 
Lamy, 2002), when the spherical symmetry is assumed to hold locally. The importance of 
the SSI method for coronal density determination has been demonstrated by wide applica- 
tions of the derived densities such as in testing models of the acceleration mechanism of 
the fast solar wind (Quemerais et al., 2007; Lallement et al., 2010), interpreting sources of 
type II and type IV radio bursts (Cho et ah, 2007; Shen et al., 2013; Ramesh et al., 2013; 
Zucca et al., 2014), and determining the coronal magnetic field strengths from fast magne- 
tosonic waves by global coronal seismology (Kwon et al., 2013). 

However, in contrast to extensive applications, few studies have evaluated the SSI 
method. Gibson et al. (1999) compared the white-light densities to those determined from 
the density-sensitive EUV line ratios of Si IX 350/342 A observed by SOHO/CDS, and 
found that densities determined from these two different analysis techniques match ex- 
tremely well in the low corona for a very symmetric solar-minimum streamer structure. 
Similarly, Lee et al. (2008) compared densities of various coronal structures determined 
by inverting MLSO MK4 pB maps and from the line ratios of Ovi 1032/1037.6 A ob- 
served by SOHO/UVCS, and found that the mean densities in a streamer derived by the 
two methods are consistent, while the coronal densities for a coronal hole and an ac- 
tive region are within a factor of two. These results are encouraging, and suggest that 
the 2D white-light density distribution in coronal structures can be very useful for other 
studies. Still, a detailed assessment is required for its better application, such as infor- 
mation about the limitation of the SSI method and the uncertainty of derived densities. 
To achieve this goal, one may use synthetic pB images from 3D densities of the corona 
reconstructed by tomographic techniques (Frazin and Janzen, 2002; Frazin et al, 2007, 
2010; Vasquez et al., 2008; Kramar et al., 2009; Butala et al, 2010; Barbey, Guen- 
nou, and Auchere, 2013), or simulated by global 3D MHD models (Mikic et al., 1999; 
Linker et al., 1999; Usmanov et al., 2000; Groth et al., 2000; Hayashi, 2005; Riley et al., 
2006; Feng, Zhou, and Wu, 2007; Hu et al., 2008; Lionello, Linker, and Mikic, 2009; 
van der Holst et al., 2014). 

The tomographic technique is a sophisticated method that reconstructs optically thin 3D 
coronal density structures using observations from multiple viewing directions. The use of 
this method in solar physics was previously proposed by Davila (1994), and later this method 
has been applied to the SOHO/LASCO (Frazin and Janzen, 2002) and STEREO/COR 1 data 
(Kramar et al., 2009, 2014). For a solar coronal tomography based on observations made 
by a single spacecraft or only from the Earth-based coronagraph, data typically need to be 
gathered over a period of half a solar rotation, so, generally, only structures that are sta- 
tionary over about two weeks can be reliably reconstructed. Therefore, this technique is not 
applicable to eclipses or, perhaps, periods of high level of solar activity, although the 3D 
coronal electron density can be routinely computed (Butala, Frazin, and Kamalabadi, 2005; 
Kramar et al., 2009). Similarly, as global MHD models of the corona need measurements 
of photospheric magnetic-field data over a solar rotation, it is also difficult to use the MHD 
method to reconstruct dynamic or rapidly evolving coronal structures that match the obser- 
vations. Thus, the SSI analysis could be very useful in some cases when the tomography 
is unsuitable. In addition, the SSI method is useful for investigating the coronal density 
variability over a long-term period (several solar cycles) when modern quality synoptic ob- 
servations were not available and relate it to the modern state-of-the-art reconstructions. 

In this study, we choose the 3D coronal density obtained by tomography as a model to 
estimate uncertainties of the SSI method. Vasquez et al. (2008) compared the tomographic 
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reconstruction and a 3D MHD model of the corona, and found that at lower heights the MHD 
models agree better with the tomographic densities in the region below 3.5 R Q , but become 
more problematic at larger heights. They also showed that the tomographic reconstruction 
has more smaller-scale structures within the streamer belt than the model can reproduce. 
Moreover, the tomographic reconstruction is entirely based on coronal observations, while 
the MHD models are primary based on photospheric boundary conditions. This suggests 
that the tomographic reconstructions are more realistic and thus may be more suitable to be 
used as a model to evaluate uncertainties of the SSI method. 

This article is organized as follows: Section 2 describes two SSI methods and their re- 
lationship. Section 3 presents the evaluation of the SSI method. We demonstrate the 3D 
density reconstruction by the SSI method based on real data in Section 4. The discussion 
and conclusions as well as potential extension of our work are given in Section 5. 


2. Spherically Symmetric Inversion Method 


Following derivations in Billings (1966), pB at a point, P, on the plane of sky (POS) can be 
expressed as 

P B{p) = — f [(l-u)A(r) + uB(r)]P-N(s)ds, (1) 

1 - u/3 \ 16 / Jlos r 

where the POS is defined as a plane crossing the Sun’s center and perpendicular to the 
LOS, N is the electron density, p is the perpendicular distance between the LOS and Sun 
center, r is the radial distance from the Sun center, and s is the distance measured from P 
along the LOS. These distances are related by r 2 = p 2 + s 2 . B 0 is the mean solar brightness, 
ay- = | Trrg = 6.65 x 1CL 25 cm 2 , is the Thomson scattering cross section for a single electron, 
where r 0 is the classical electron radius. Note that r, 2 = 7.94 x 10~ 26 cm 2 is the Thomson 
scattering cross section as referred to in Billings (1966). u is the limb-darkening coefficient. 
In addition, A(r) and B(r ) are geometrical factors given by Billings (1966), 


A(r) = cos £2 sin” £2 , 




, , < 1 + 3 sin £2 

1 — 3 sin" £2 — cos" £2 , 

sin £2 


In 


1 + sin £2 
cos £2 


(2) 

(3) 


where the angle £2 is defined by sin £2 = R Q /r. 

If the electron density is a function of r only, Equation (1) can be written in the form 


pB{p) = C 



[( l-u)A(r) + uB(r)]N(r ) 


p 2 dr 

rs/r 2 - p 1 


(4) 


where r and p are in units of B 0 , and C = (3/8)ct7’/? 0 B 0 /(1 — m/ 3). Van de Hulst (1950) 
developed a technique for inverting pB measurements (called the Van de Hulst inversion). To 
implement this technique, one needs first to fit the pB data using a curve in the polynomial 
form, specifically. 


pB(p) = ^a k p k . 

k 


(5) 


After the coefficients a k are determined, the solution of Equation (4) is given by 

N(r) = EM/c k+l )r~^ 

C(S\-u)A{r) + uB{r)Y 


where 


(6) 
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Ck ■ 


r n/2 

Jo 


*... V5Fr((Jfc+l)/2) 

COS (p cicp = — . 


(7) 


F(k/2 + 1) 

where T is the gamma function. Note that the constant C here is different from that in Van de 


Hulst (1950), because we have used the expression of pB in Billings (1966). More generally, 
the index k of radial power law in Equation (5) is real but not necessarily integer valued. This 
leads to the modified Van de Hulst technique assuming pB(p) = ^2 k a k p~ dk , where a k and 
d k are the fit coefficients, then Ck + i and r _(A:+1) in Equation (6) need to be replaced with ca k+ 1 
and r~ < - dk+l ' 1 ( e.g ., Saito, Poland, and Munro, 1977; Guhathakurta, Holzer, and MacQueen, 
1996; Gibson et al., 1999). 

Hayes, Vourlidas, and Howard (2001) developed another SSI technique by assuming the 
radial electron density distribution in the polynomial form, N(r ) = bkr~ k , and have used 
it to invert total brightness observations. We here adopt this technique for the pB inversion 
and call it the spherically symmetric polynomial approximation (SSPA) method for short. 
By substituting the polynomial for N(r) in Equation (4), we obtain 


P B{p) = Y J bkG k {p ), 

k 


( 8 ) 


where 

G k (p) = C [(1 - u)A(r) + uB(r)](r- k )^J- . (9) 

J P r^/r 2 - p- 

For easier calculation, using the substitution r = p/cos# the integral can be transformed 
into 


|.ir/2 

G k (p) = Cp~ k+i [(l-u)A(r(6)) + uB(r(0))]cos k 6d6, (10) 

Jo 

where 9 is the angle from the POS to a direction from the Sun center to a point of distance s 
along the LOS. Since the integral G k {p) can be numerically calculated for all desired impact 
distances p and exponents k, substituting the observed pB curve along a radial trace for the 
left-hand side of Equation (8) the only unknowns are the coefficients b k . We determine the 
coefficients by a multivariate least-squares fit to the curve of pB(p) (using the function of 
svdfit provided by IDL, Interactive Data Language). The radial distribution of the electron 
density is then obtained by substituting the resulting coefficients directly into the polynomial 
form. To select appropriate degrees of polynomial we performed some experiments using 
the coronal density models for the region between 1 .5 and 6 R,., in Saito, Poland, and Munro 
(1977), and found that choosing the first five terms (k = 1-5) can determine N (r) with the 
relative errors within 5 % and reproduce pB measurements with the relative errors within 
1 %. We therefore use the five-degree polynomial fits for all the SSPA inversions 1 in our 
study. 

To determine the relationship between the SSPA and Van de Hulst inversions theoreti- 
cally, we use Taylor series to approximate the functions A{r) and B(r) in the case when 
r > 1 with r in unit of R e , 


A(r ) 
B(r ) 



2 / 1\ 2 4 / 1 \ 4 3 / 1 \ 6 

3 W ~ 15 W “ 40 V/ + 


( 11 ) 

( 12 ) 


*IDL codes of the SSPA method for inversions of COR1 pB data in ID and 2D are available for download 
(http://solar.physics.montana.edu/wangtj/sspa.tar). 
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For the Taylor polynomial approximations of degree six above, the relative errors are lower 
than 1 % when r > 1.5. Since A(r) and B(r ) are on the order of 0(l/r 2 ) for very large 
r, by keeping only the terms of 1 / r 1 we have (1 — u)A(r) + uB{r) (1 — u/3)(l/r 2 ), 
which corresponds to the point source approximation (Frazin et al., 2010). Applying this 
approximation to Equation (6) for the Van de Hulst inversion, we obtain 


N(r) = 


C( 1 



(13) 


This verifies the polynomial form of N (r) assumed in the SSPA method. Conversely, if N (r) 
is given in the same form as Equation (13), likewise we can recover pB(p) = '}2 k akP~ k as 
assumed in the Van de Hulst inversion from the SSPA inversion using Equations (8) and (10) 
with the point-source approximation for (1 — u)A(r) + uB(r). Thus, we have proved that 
the SSPA and Van de Hulst inversions are identical on the order of 0(l/r 2 ), and their 
difference in higher orders can be reduced by increasing the degrees of k. Note that the 
above analysis is also valid for the case when k is real. In Section 3.2, we show that the 
coronal densities determined from COR1 images by the SSPA and Van de Hulst inversions 
with k = 5 are almost identical. Therefore, the assessment of the SSPA inversion in the 
following is equivalent to that of the Van de Hulst inversion. 


3. Evaluation of the SSI Method 

3.1. Coronal Density Model for Evaluating the SSI Method 

As a model for evaluating the SSI method we used a 3D coronal electron density obtained 
by tomography method applied to white-light coronagraph data (Kramar et al, 2009). The 
data were acquired by the inner coronagraph (COR1) telescopes onboard the twin Solar 
Terrestrial Relation Observatory (STEREO) spacecraft. COR1 observes the white-light K 
corona from about 1.4 to 4 R Q in a waveband of 22.5 nm wide centered on Ha line at 656 
nm with a time cadence of five minutes. The regularized tomographic inversion method by 
Kramar et al. (2009) with the limb-darkening coefficient u = 0.6 provided the reconstruction 
of a 3D coronal electron density for the period of 1 - 14 February 2008 (consisting of 28 pB 
images from COR1-B) that corresponds to Carrington Rotation (CR) 2066. The scattered 
light in the pB data was removed by subtracting a combination of the monthly minimum 
and the roll minimum backgrounds (Thompson et al., 2010). The reconstruction domain is 
a 128 3 rectangular grid covering a spherical region between 1.5 and 4 R Q . We evaluate the 
SSPA inversion method using this 3D density reconstruction as a model. 

To synthesize pB images observed by COR1-A and -B at a given time, we first deter- 
mined the orbital positions of the spacecraft in Carrington coordinates using the routine 
get_stereoJonlat provided by SolarSoftWare (SSW). The 3D density grid was then trans- 
formed from the Carrington heliographic system to the projected coordinate system viewed 
from STEREO-A or -B (Thompson, 2006). Therefore, the LOS integral of N (r) in Equation 
(1) becomes a simple summation in the ^-direction (defined along the LOS). Figure 1 shows 
the tomographically reconstructed 3D coronal electron density at 2 R Q . Figure 2 shows the 
density distributions of the corona in the POS for COR1-A and -B at 12:00 UT on 8 Febru- 
ary 2008. Figure 3 shows the corresponding synthetic pB images, which represent the ideal 
measurements without contamination by the scattered light and instrumental noises. 
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Ne for CR 2066 of 2.0 R 0 



0 60 120 180 240 300 10* cm"’ 

Carrington longitude (deg) 


Figure 1 Spherical cross section of the tomographically reconstructed 3D coronal electron density for CR 
2066 at a heliocentric distance of 2.0 Rq. Two vertical solid lines marked with AE and AW indicate the 
positions of the east and west limbs in the COR1-A image observed at 12:00 UT on 8 February 2008, and the 
dashed lines with BE and BW denote the positions of east and west limbs in the COR1-B image at the same 
time. The diamonds and pluses indicate the locations for the SSPA inversions shown in Figures 7 and 8. 


(a) Tomographic density for C0R1— A (b) Tomographic density for C0R1-B 
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Figure 2 Cross sections of the tomographic 3D coronal density in the plane of the sky (POS), corresponding 
to positions of COR1-A (a) and COR I B (b) at 12:00 UT on 8 February 2008. 


3.2. Radial Distribution of Electron Densities in Solar-Minimum Streamer 

Coronal streamers are the most conspicuous, large-scale structures in the extended corona. 
The streamer belt is a usually continuous sheet of enhanced density associated with the mag- 
netic neutral sheet or the current sheet (Schulz, 1973; Guhathakurta, Holzer, and MacQueen, 
1996; Saez et al., 2007; Kramar et al., 2009, 201 1). Its shape is progressively deformed from 
a rather flat plane at minimum solar activity to a highly warped surface at maximum solar 
activity ( e.g ., Hu et al., 2008). Some previous studies have suggested that the SSI assump- 
tions are suitable to symmetric streamers at low latitude, in particular, the solar-minimum 
streamer belt (e.g., Guhathakurta and Fisher, 1995; Gibson and Bagenal, 1995; Gibson et al., 
1999). Here we first test the SSPA inversion of the solar minimum streamer belt. A streamer 
belt during CR 2066 corresponding to solar minimum of Solar Cycles 23/24 is shown in 
Figure 1. 

For instance, we used the synthetic pB images at 12:00 UT on 8 February 2008 when 
the separation angle between STEREO- A and -B was about 45°. Figure 3 shows that the 
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Figure 3 Synthetic pB images based on the 3D coronal density model reconstructed by tomography, 
(a) Viewed from COR1-A, and (b) from COR1-B at 12:00 UT on 8 February 2008. Carrington longitudes 
of the viewing direction are marked at the bottom of the images. Solid lines mark the positions where the 
pB intensity profiles shown in Figure 4 are extracted. Short bars with a size of 0.3 1 Rq show the scale over 
which the pB intensity is averaged across the radial cut. Four circles (dotted lines) at 1.6, 2.0, 2.5, and 3.0 Rq 
mark the paths along which the density profiles are shown in Figure 1 1 . 


Figure 4 Radial profiles of pB 
for COR1-A (solid line) and 
COR 1 -B (dashed line), extracted 
along the radial cuts shown in 
Figure 3. The error bars are the 
standard deviations for an 
average over five pixels across 
the cut. 


1.5 2.0 2.5 3.0 3.5 4.0 

Radio distance (R @ ) 



streamer belt on the east limb is almost edge-on, as inferred from the similar appearance in 
COR1-A and -B. In the edge-on condition, the coronagraph is looking along the streamer 
belt; i.e., where all the streamers are at the same latitude behind each other along the LOS. 
In contrast, the streamer belt on the west limb in COR1-B is seen face-on, showing a dis- 
tinctly different shape from that in COR1-A. We chose a radial trace near the middle of the 
east-limb streamer (see marked positions in Figure 3), where the pB profiles for COR1-A 
and -B are nearly identical (Figure 4), and assumed that this location best meets the SSI con- 
dition. We then derived the radial distribution of the electron density by fitting the pB data 
along this radial trace between 1.6 and 3.9 Rq using both the SSPA method and the Van de 
Hulst technique (pbjnverter with k = 5 in SSW) for an illustrated comparison. We show the 
inversion results in Figure 5. Note that since COR1 does not perform well at p < 1.55 Rq at 
some position angles because of interference from the occulter (see Frazin et al., 2012), this 
led to defective pB signals and thus unreliable reconstructions within this region, therefore 
we constrain all SSPA inversions in the following to regions with p > 1.6 Rq. 

The electron densities obtained from both pB inversion methods agree well in radial 
distribution with the model 3D densities in the POS for selected longitudes corresponding 
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(a) pB inversion for C0R1— A 




(c) pB inversion for C0R1-A (d) pB inversion for C0R1-B 



Figure 5 (a) and (b): The synthetic pB data (denoted with pluses) along a radial cut shown in COR1-A and 
-B images in Figure 3, and their least-squares fit (solid line) using the SSPA method, (c) and (d): The density 
profiles derived with the SSPA inversion (thin solid line) and the Van de Hulst inversion (dashed line), as well 
as the model densities in the POS (thick solid line). Note that in panel (d) the density profiles obtained by the 
SSPA and the Van de Hulst inversions are overlaid completely. For comparison, the pB profile and the density 
model for the equatorial background derived by Saito, Poland, and Munro (1977) from the Skylab data are 
plotted with dot-dashed lines. In (c) and (d) the overplotted curves also include the solar-minimum streamer 
densities (dotted line) obtained by Gibson et al. (1999) from Mauna Loa/LASCO C2 data. 


to the edge-on streamer positions (see Figures 5(c) and (d)). This confirms our theoretical 
prediction that the PPSA and the Van de Hulst inversions are equivalent. The average ratio 
of the SSPA to the model densities along the selected radial cut between 1.6 and 3.0 Rq is 
0.86 ± 0.06 and 0.95 ± 0.10 for COR1-A and -B, respectively. Note that the weak increase 
in density with height between 3.5 and 4.0 R Q may be due to FOV effects (see discussion 
in Section 5). In addition, the lower panels of Figure 5 show that both the tomographic and 
SSPA density profiles determined for streamers are similar (with differences of ~ 20 %) 
to those previously obtained by SSI method in the solar minimum (Gibson et al., 1999; 
Saito, Poland, and Munro, 1977). 

To examine to what extent the SSPA solution is consistent with the 3D density model, we 
compared the SSPA density profile with the model density profiles in 13 angular sections 
with the POS within 30° in Figures 6(a) and (b), which shows that the SSPA solution is clos- 
est to the distribution of 3D densities in the POS. This is expected because the pB integrals 
along the LOS are most heavily weighted toward the regions near the POS (see Quemerais 
and Lamy, 2002; Frazin et al., 2010). 
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Figure 6 Comparison of the SSPA density profile (dashed line) with the model densities in 13 angular 
sections (solid lines). The thick black line stands for the case in the POS, while thin color lines plot cases with 
angles relative to the POS to be ± 5, 10, 15, 20, 25, and 30 degrees, where the negative value represents an 
angle measured inward and the positive value an angle measured outward, (a) For COR1-A. (b) for COR1-B. 
/,,\F and !-]][■ marked in (a) and (b) represent the longitude of the east limb in COR1-A and -B. 


3.3. Longitudinal and Latitudinal Dependence of the Streamer Density 

We evaluate the SSPA inversion of the streamer as a function of longitude in the following. 
This corresponds to making the SSPA inversion of synthetic pB images at different “virtual 
observing times”. We synthesized the pB images at nine times for COR1-A and -B in the 
period from 20:00 UT, 6 February to 04:00 UT, 10 February with intervals of 10 hours, 
centered at 12:00 UT on 8 February 2008 (the time for the instance analyzed in the last 
section). The Sun rotated by 44° during this period, approximately equal to the separation 
angle between STEREO-A and B. So the ending-time location analyzed in COR1-A is ap- 
proximately superposed with the starting-time location in COR1-B as shown in Figure 1. 
We fit the pB data along the same radial trace between 1.6 and 3.9 Rq at these nine times 
using the SSPA method. Figures 7 and 8 show comparisons of the obtained density pro- 
files with the model densities in the POS for COR1-A and COR1-B, respectively. For the 
analyzed locations with longitude between 60° and 106° (approximately between the limbs 
BE and AE as shown in Figure 1), the density profiles derived by the SSPA method have a 
similar shape as the model profiles over almost the whole FOV range (1.6 -3. 8 R 0 ), but the 
magnitudes are smaller than those of the model by about 20 %-50 %. The reason for the 
good agreement during this period might be that the streamers in the LOS at the analyzed 
location are near the POS. The model density increases near the edge of the FOV seen in 
Figures 7(c) -(e) and Figures 8(g) — (i) most likely result from the FOV effects in the to- 
mographic inversion (Frazin et al., 2010). Figure 9 shows the SSPA density as a function 
of longitude at five heights (1.6, 2.0, 2.5, 3.0 and 3.5 R Q ). The comparison with the model 
confirms the result above that the better agreement lies at the locations with the longitude 
in 60° - 106°. Moreover, the agreement tends to be better at lower heights. The average dif- 
ference between the SSPA and model densities are ~ 20 % for the region between 1.6 and 
2.0 Rq, while ~ 40 % - 50 % for 2.5 - 3.5 Rq . The reason might be that the SSI condition is 
better met at lower heights where the streamer has a larger extent in the LOS, and thus the 
LOS integration in Equation (10) over a very long distance (compared to the width of the 
streamer) becomes more reasonable. 
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Figure 7 The electron density profiles derived by the SSPA inversion (dashed line) as a function of longitude. 
The analyzed radial trace is located at the position with Carrington latitude of 9° and longitude marked at the 
top of each panel and shown with pluses in Figure 1 . The solid line in each panel represents the corresponding 
model densities in the POS. 


Now we evaluate the latitudinal dependence of SSPA inversions of the corona. By as- 
suming that the SSI condition holds locally for all angular positions around the Sun, the 2D 
coronal density can be derived from a pB image by fitting the radial profile at each angular 
position using the SSPA method. For the case analyzed in Section 3.2, we reconstructed 
2D coronal densities in the POS by fitting radial profiles between 1.6 and 3.7 R Q for 360 
position angles with intervals of 1° from both synthetic and observed pB images, and show 
the results in Figure 10. Compared with the tomographic densities in the POS (see Fig- 
ure 2), the SSPA coronal densities determined from the observed images do not show very 
low-density regions around streamers, while those from the synthetic images have a small 
region of zero (or negative) densities at the southern side of the west-limb streamer, where 
the synthetic pB signals are very low (below the background level in non-streamer regions). 
The zero-density regions in the tomographic reconstructions could be caused by either (or 
both) coronal dynamics or a real very low density in these coronal regions (see discussion 
in Section 5). 
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Figure 8 Same as Figure 7, but for C0R1-B. Locations of the analyzed radial trace are marked with dia- 
monds in Figure 1, whose latitude is 9°, and the longitude is marked at the top of each panel. The solid line 
represents the model densities in the POS, the dashed line represents the SSPA densities. 




Figure 9 The coronal densities at 1.6, 2.0, 2.5, 3.0, and 3.5 Rq as a function of Carrington longitudes, 
derived with the SSPA inversion (dashed lines) in comparison with those by the model in the POS (solid 
lines). The Carrington latitude of the analyzed radial trace is 9°. (a) For COR1-A. (b) For COR1-B. The 
vertical dotted lines indicate the middle position of the analyzed longitude range. 
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(a) SSPA density from synthetic pB for C0R1 -A 
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Figure 10 The 2D coronal density derived using the SSPA method from the synthetic pB images, corre- 
sponding to positions of COR1-A (a) and COR1-B (b) at 12:00 UT on 8 February 2008. Carrington longi- 
tudes of the viewing direction are marked at the bottom of the density maps, (c) and (d): same as (a) and (b), 
but derived from the real observed pB images at the same time. 


For a quantitative comparison, we plot the SSPA and model density profiles as a function 
of position angles at four heights (1.6, 2.0, 2.5, and 3.0 R Q ) in Figure 11. We measured 
the peak densities and angular widths (in FWHM) of the three edge-on streamers marked 
SI - S3. The ratio of peak densities from the SSPA to the model is on average 0.82 ± 0.06, 
0.72 ± 0.09, and 0.62 ± 0.12, and the ratio of angular widths is on average 1.24 ± 0.02, 
1.67 ± 0.12, and 1.90 ± 0.27 at 2.0, 2.5, and 3.0 Rq, respectively. The results indicate that 
the SSPA inversion underestimates the peak density of streamers by about 20 %-40 % 
and overestimates the angular width by about 20 % -90 %. The increase of deviations with 
height in the peak density of streamers by the SSPA method from the model may be caused 
by the narrowing of the streamer (belt) with height (see Figure 2 of Kramar et al., 2009). 
As the streamer width along the LOS decreases with the height, the condition for the SSI 
assumption becomes worse. The reason for the latitudinal spreading of the SSPA density rel- 
ative to the model may be that the analyzed streamers are seen not exactly edge-on, i.e., the 
streamer belt at the analyzed locations is actually tilted somewhat away from the equatorial 
direction (or the LOS near the limb, see Figure 1). In this case at the locations near the edge 
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Figure 11 The coronal densities along four circular paths at 1.6 (a), 2.0 (b), 2.5 (c), and 3.0 Rq (d) for 
COR1-A. The position angle is counted anticlockwise from the north pole. The eastern part of each circle 
(position angles 0°-180°) has a longitude of 106.2°, its western part a longitude of 286.2°. The solid line 
represents the model densities, the dashed line and the dotted line represent the SSPA densities from the 
synthetic and observed pB images, respectively, (e)-(h): same as panels (a)-(d), but for COR1-B. The 
eastern part of each circle has a longitude of 60.5°, its western part a longitude of 240.5°. 


of streamer in the pB image the contribution is larger from points along the LOS behind or 
in front of the POS, leading to an overestimation of the model density in the POS by the 
SSPA inversion at these places. In addition, we note for the face-on streamer (marked S4) in 
COR1-B that its SSPA density profiles are consistent with the model densities as well. This 
is because this streamer is by coincidence located near the POS at this instance (Figure 1). 

In Figure 1 1, we also compare the SSPA densities inverted from the synthetic pB image 
with those from the observed pB image. They are generally consistent, except in the region 
near the occulter of COR 1 , where the SSPA inversions from the observed pB are much larger 
(~ 2 - 3 times) than the model densities (see panels (a) and (e) in Figure 11). The version of 
tomographic reconstruction used here most likely underestimated the density there because 
the solution at the grid points close to the occulter is less constrained by the observational 
data. 


4. Reconstruction of the 3D Coronal Density Using the SSPA Method for CR 2066 

In the above sections, we evaluated the SSPA method by comparing the derived density dis- 
tribution as a function of radial distance, longitude, and latitude with the 3D tomographic 
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model. The assessments indicate that the SSPA inversion can determine the 2D coronal den- 
sity from a pB image, which approximately agrees with the 3D tomographic densities for 
streamers near the POS. This suggests that we may reconstruct a 3D density of the corona 
by applying the SSPA inversion to a 14-day data set from COR1-A or -B. In this section, 
we demonstrate the SSPA 3D density reconstructions using the same data set (consisting of 
28 pB images from COR1-B) as used for the reconstruction of 3D tomographic model, and 
the simultaneous COR1-A data set. The data cadence of about two images per day corre- 
sponds to a longitudinal step of about 6°. First we determined the 2D density distribution by 
fitting the radial pB data between 1.6 and 3.7 R Q using the SSPA inversion at 142 angular 
positions (with intervals of 2.5°) surrounding the Sun for each image. Then we mapped the 
east-limb and west-limb density profiles of all images to make a synoptic map at a certain 
height, based on their Carrington coordinates (neglecting the inclination of the solar rota- 
tional axis to the ecliptic). A 3D density reconstruction was made of 25 synoptic maps for 
the radial heights from 1.5 to 3.9 Rq with an interval of 0.1 R Q . For each synoptic map, we 
converted the irregular grid into the regular grid using the IDL function, trigrid. Although the 
3D density reconstruction can be made using the SSPA method from the COR1 data with 
a cadence as high as five minutes, the higher temporal resolution may not help to improve 
its actual angular resolution in longitude because of the intrinsic limitations of the SSPA 
method (see Appendix). 

Figure 12 shows the SSPA-reconstructed 3D coronal density at 1.6, 2.0 and 2.5 Rq for 
COR1-A and -B. Some discontinuities of the density are seen at two longitudes that separate 
the two regions made of the east-limb and west-limb inversions. These flawed structures may 
be due to temporal changes of streamers and/or the effect due to neglecting the solar axial 
tilt in the reconstruction. We smoothed these discontinuities by averaging the two recon- 
structions from COR1-A and COR1-B and then made a 10° x 10° smoothing. In Figure 13, 
we show the ratio of tomographic density to the SSPA average density for the COR1-A and 
-B reconstructions with smoothing for two cases, one for the SSPA densities obtained using 
the synthetic data (panels (a) -(c)), and the other for the SSPA densities obtained using the 
real data (panels (d)-(f)). Both cases indicate that the density ratios in the streamer belt 
are very close to 1, within of about a factor of two (~0.5-l at 1.6 Rq, ~ 1-2 at 2.0 R Q , 
and ~ 1 - 3 at 2.5 Rq). The good consistency validates that the SSI assumptions are very 
appropriate for the streamer belt in the solar minimum. 

For a quantitative comparison between the tomography and the SSPA densities (obtained 
from real data), we also show their density profiles along an equatorial cut in Figure 14. 
They match best at 2 R Q . The same result is indicated by the pixel-to-pixel scatter plots 
in Figure 15. We obtain a ratio of the SSPA to the tomographic average density of 1.02, 
and their linear Pearson correlation coefficient is 0.93 for the reconstructions at 2 R Q . In 
addition, the scatter plots show that the dispersion is larger for those pixels with lower values, 
where the SSPA densities are much higher than those obtained by the tomography. 


5. Discussion and Conclusions 

In this study, we have, for the first time, evaluated the SSI method using the 3D model of 
the coronal electron densities reconstructed by tomography from the STEREO/COR 1 pB 
data. Our study is instructive for a more efficient use of the SSI technique to invert the pB 
observations from ground- and space-based coronagraphs, in particular, the COR1 data. We 
demonstrated both theoretically and observationally that the SSPA method and the Van de 
Hulst inversion are equivalent SSI techniques when the radial densities or pB signals are 
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(a) The SSPA density at 1.6 R 0 for C0R1-A 
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Figure 12 The SSPA-reconstructed 3D coronal densities from the two-week dataset of COR1 pB images. 
(a)-(c): spherical cross sections at 1.6, 2.0, and 2.5 Rq for COR1-A, and (d)-(f) for COR1-B. In panels 
(a), (b), (d), and (e), the dotted and the dashed line represent the positions of the east and the west limb, 
respectively, of the earliest-observed pB image (at 00:25 UT, 1 February 2008). In panels (c) and (f), the 
vertical dotted and dashed lines represent the positions of the east and west limb, respectively, for all pB 
images used in the reconstructions. 


assumed in the polynomial form of high degrees (more than two). The polynomial degree of 
five is suitable for COR1 data inversions. Thus, assessment results of the SSPA method can 
also be applied to the Van de Hulst inversion technique. We determined radial profiles of the 
streamer density from the COR1-A and -B synthetic pB images as well as their longitudinal 
and latitudinal dependencies. We found that the SSPA density values are close to the model 
for the core of streamers near the POS, with differences (or uncertainties if we regard the 
model input as a true solution) typically within a factor of ~ 2. This result is consistent with 
those evaluated using UV spectroscopy (Gibson et al, 1999; Lee et al., 2008). We found 
that the SSPA density profiles tend to better match the model at lower heights (< 2.5 R Q ). 
Our results confirm the suggestion made in some previous studies that the SSI assumption is 
appropriate for the edge-on streamers or the streamer belt during solar minimum. We suggest 
that the edge-on condition for streamers may be determined by the tomography method or 
by examining the consistency between simultaneous pB measurements from COR1-A and 
-B in radial distribution, when the two spacecraft have a small angular separation (e.g., less 
than 45°). We also found the SSPA streamer densities extend more in both longitudinal and 
latitudinal directions than in the model. 

We demonstrated the application of the SSPA inversion for reconstructions of the 3D 
coronal density near solar minimum, and showed that the SSPA 3D density for the streamer 
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Figure 13 The ratio of tomographic density to the SSPA average density for COR1-A and -B at 1 .6, 2.0, and 
2.5 Rq for two cases, (a) -(c): for the SSPA 3D coronal density obtained using synthetic pB data based on 
the 3D density model, (d)-(f): for the SSPA 3D coronal density obtained using the real pB data observed by 
COR1-A and -B. The colors represent the ratio values in a range 0-5 (see color bars at the bottom). A version 
with continuous color scale is available as electronic supplementary material. 


belt is roughly consistent in both position and magnitude with the tomographic reconstruc- 
tion. The synoptic density maps derived by the SSPA method show some discontinuities at 
the longitudes that separate the regions made of the east- and west-limb inversions. These 
discontinuities may be due to temporal changes of coronal structures and/or the effect of 
tilt of the solar rotation axis on the poles’ visibility, which was not considered. They can be 
smoothed out during post-processing by smoothing and averaging the density distributions 
from COR1-A and -B, but such treatments will reduce the spatial resolution. In comparison, 
the tomographic inversion can fully take into account the tilt effect of the solar pole and pro- 
duce smoother density distributions in these discontinuity regions (Frazin and Janzen, 2002; 
Kramar et al., 2009). We estimate that the SSPA method may resolve the coronal density 
structure near the POS with an angular resolution of ~ 50° in longitudinal direction. Given 
this limitation, the SSPA reconstruction using pB data with higher cadence ( e.g ., more than 
three images per day) would not help to improve its actual angular resolution. 

Although the current state of the tomographic method has allowed us to routinely obtain 
the 3D coronal densities, we speculate that the SSI method could be complementary to the 
tomography when used for the interpretation of observations for instance during of solar- 
activity maxima, in some regions where tomography gives zero-density values, or in the 
regions near the edges of the FOV. 
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Figure 14 Comparison of the 
SSPA density profiles with the 
tomographic density profiles 
(thick solid line) along the 
equator at (a) 1 .6, (b) 2.0, and 
(c) 2.5 Rq. The dotted and the 
dashed lines represent the SSPA 
coronal densities for COR1-A 
and -B, respectively, while the 
thin solid line represents their 
average with a 10-degree 
smoothing. 
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Figure 15 The pixel-to-pixel scatter plots of the SSPA 3D coronal density (averaged for COR1-A and -B 
with a 10-degree smoothing) versus the tomographic 3D coronal density, at (a) 1.6, (b) 2.0, and (c) 2.5 Rq. 
The linear Pearson correlation coefficient (px,y) and the ratio of the average densities ((Y) / (X)) are marked 
in each plot. 


The zero-density values in the tomographic reconstruction (so-called zero-density arti- 
facts, ZD As) might be caused by either (or both) coronal dynamics (Frazin and Janzen, 2002; 
Frazin et al., 2007; Vasquez et al., 2008; Butala et al., 2010) or (and) real very low density 
in these coronal regions, which are below the error limit in the tomography method (Kramar 
et al., 2009). The latter reason is also supported by results of MHD modeling (Airapetian 
et al., 2011). In the former case, the SSI method could be complementary to tomography, 
while in the latter case the SSPA method gives much higher values than in the tomographic 
model. The FOV artifacts in the tomography are due to the finite coronagraph FOV that 
causes the reconstructed density to increase at the regions close to the outer reconstruction 
domain (Frazin et al., 2010; Kramar et al., 2009). However, these FOV artifacts can be 
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reduced by extending the outer reconstruction domain beyond the coronagraph FOV limit 
(Frazin et al., 2010; Kramar et al., 2013). Another way to obtain more correct density values 
in this region might be by using the SSPA method, which does not imply strict outer bound- 
aries for LOS integration. Thus the estimate of uncertainties of the SSPA method should be 
limited to the regions with distances smaller than about 3.5 Rq for the tomographic model 
used. An MHD model to produce artificial data can be useful for this test as well. However, 
this will be a subject of future research. 

Tomography generally assumes that the structure of the corona is stable during the 
observational interval, e.g., two weeks of observations made by a single spacecraft, al- 
though for some coronal regions that are exposed to the spacecraft for only about a week 
during the observation the stationary assumption can be reduced to about a week (Kra- 
mar et al., 2011). However, this assumption is hard to meet during solar maximum or 
times of enhanced coronal activity. The CME catalog in the NASA CDAW data cen- 
ter shows that the CME occurrence rate increased from ~ 0.5 per day near solar min- 
imum to ~6 near solar maximum during Solar Cycle 23 (Gopalswamy et al., 2003; 
Yashiro et al., 2004). Although our assessment results for the SSPA method are based on a 
static coronal model, their validity may not be limited to the static assumption, because the 
key factor for the SSPA inversion for obtaining a good estimate of the 2D coronal density is 
an instantaneous (local) symmetric condition for coronal densities along the LOS. The min- 
imal size of this local symmetry is limited by the angular resolution in longitude, which is 
about 50°. Therefore, the SSI method can be used to estimate the density of a dynamic coro- 
nal structure in terms of weighted average over the region with this angular size in longitude. 
For this reason it may be a better choice to use a combination of tomography and the SSI 
inversions for interpretations of radio bursts and shocks produced by CMEs in the case when 
coronal structures of interest evolve quickly with time (Lee et al., 2008; Shen et al., 2013; 
Ramesh et al., 2013). In addition, the SSI method is also often applied to the cases when 
observational data are unsuitable for tomography, e.g., solar eclipses. 

The 3D MHD models of the corona using the synoptic photospheric magnetic held data 
have been successfully used to interpret solar observations, including total eclipses and 
ground-based {e.g. MLSO/MK4) images of the corona (e.g.. Linker et al., 1999; Mikic etal., 
1999). We suggest that the evaluation of the SSI method based on such a global MHD model 
may be necessary in the future. The profits using such MHD models to estimate the uncer- 
tainties of the SSI method could be in avoidance of ZDA and FOV effects. The modeled 
corona also allows evaluations of the SSI method down to very low heights (e.g., the region 
between 1.1 and 1.5 Rq), where the corona is much more structured. Moreover, a simulated 
time-dependent corona from a time-evolving MHD model would allow us to estimate the 
uncertainties in tomography and the SSI inversion when they are applied to the dynamic 
corona, especially during solar maximum. This needs detailed investigations in the future. 
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Appendix: Estimates of Angular Resolution of the SSPA Method in Longitude 

To determine the angular resolution of the SSPA method in the longitudinal direction, a 
numerical experiment was performed using 2D coronal models. We constructed a 2D density 
model by first using the Saito, Poland, and Munro (1977) equatorial background density 
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Figure 16 2D density models of the corona from 1 .5 to 4.0 Rq in the equatorial plane with two structures 
(marked A and B) of (a) the step profile and (b) the Gaussian profile in longitude. The horizontal dotted 
line represents the POS. (c)-(e): comparisons of the SSPA density (dashed line) with the model density 
(solid line) in the POS at 2 Rq as a function of longitude for the two structures with angular separation 
of 40° (c), 50° (d), and 60° (e). The black curves represent the model density with the step profile, the red 
curves the model density with the Gaussian profile, (f): comparison between the SSPA and model densities 
at different heights. The black, red, and green dashed lines represent the SSPA density profiles at 2, 3, and 
4 Rq, respectively. Note that the model density profiles (in solid line), normalized to the coronal background, 
are the same for the three heights, (g): structures with different widths of 15° (red curves), 25° (black curves), 
and 35° (green curves). The solid lines represent the model density in the POS, the dashed lines the SSPA 
density, (h): structures with different density contrast ratios to the background of 5 (red curves), 10 (black 
curves), and 20 (green curves). The solid lines represent the model density in the POS, the dashed lines the 
SSPA density. In (f)-(h), except for one parameter that is set as different values, the other parameters are 
set to be 50° for angular separation, 2 Rq for heliocentric height, 25° for structure width, and 10 for density 
contrast ratio. The vertical dotted lines in (c) - (h) mark the central position in the structure. 
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model to build a background corona of rotational symmetry in the equatorial plane, and 
then inserting two structures into it. The structures have an angular width of 0, a density 
contrast ratio of d to the background, and a longitudinal profile following a step function or 
Gaussian function. The FWHM of the Gaussian-type structure was set to 0. Figures 16(a) 
and (b) show the two types of 2D coronal models with about the same FOV as COR1, where 
two structures are separated by an angle of 20, thus the width of the gap between them is 
also equal to 0 in the step-profile case. We assumed that the two structures (marked A and B) 
are located at longitudes of 0 and —0, respectively, that is, we defined their middle position 
as the origin of the longitude. For the cases shown in Figures 16(a) and (b), the origin of the 
longitude is just located in the POS at the west limb. 

We synthesized pB data for the west-limb region from 1.5 to 4 R e at longitudes in the 
range from — 90° to 90° using Equation (1), and then derived the electron density using 
the SSPA method by fitting the synthetic pB data. Panels (c) - (e) show comparisons of the 
SSPA density with the model density in the POS at 2 R Q as a function of longitude for 
different angular distances between two artificial structures with a density contrast ratio 
d = 10. For both types of the structure (in the step or Gaussian profile), the smallest resolv- 
able distance (or angular resolution) is about 50°. To examine the effects of radial distance, 
structure width, and density contrast on the obtained resolution, we made a parametric study, 
and show the results in panels (f)-(h). The resolution is only slightly dependent on the ra- 
dial distance and structure width, which is better for the lower heights and narrower widths, 
but almost independent of the density contrast. In addition, we also found that the obtained 
SSPA peak density is about a half of the model density for the analyzed structures in most 
of the cases, which is consistent with our results for edge-on streamers. These numerical 
experiments suggest that the SSPA solutions are more accurate for coronal structures with 
a smoother profile and larger extension in the longitudinal direction and with lower density 
contrast to the background, which can be regarded as better conditions meeting the spheri- 
cally symmetric assumption. 
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